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ABSTRACT 



Context. Protoplanetary disks are observed to remain dust-rich for up to several million years. Theoretical modeling, 
on the other hand, raises several questions. Firstly, dust coagulation occurs so rapidly, that if the small dust grains are 
not replenished by collisional fragmentation of dust aggregates, most disks should be observed to be dust poor, which is 
not the case. Secondly, if dust aggregates grow to sizes of the order of centimeters to meters, they drift so fast inwards, 
that they are quickly lost. 

Aims. We attempt to verify if collisional fragmentation of dust aggregates is effective enough to keep disks 'dusty' by 
replenishing the population of small grains and by preventing excessive radial drift. 

Methods. With a new and sophisticated implicitly integrated coagulation and fragmentation modeling code, we solve 
the combined problem of coagulation, fragmentation, turbulent mixing and radial drift and at the same time solve for 
the 1-D viscous gas disk evolution. 

Results. Wc find that for a critical collision velocity of 1 m/s, as suggested by laboratory experiments, the fragmentation 
is so effective, that at all times the dust is in the form of relatively small particles. This means that radial drift is small 
and that large amounts of small dust particles remain present for a few million years, as observed. For a critical velocity 
of 10 m/s, we find that particles grow about two orders of magnitude larger, which leads again to significant dust loss 
since larger particles are more strongly affected by radial drift. 
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1. Introduction 

The quest for a comprehensive understanding of the for- 
mation of planets in the dusty circumstellar disks around 
many young stars has gained significant impetus. This is 
partly because of developments in the numerical modeling 
of these processes, but also because of the high-quality in- 
frared and (sub-)millimeter data that are now available of 
hundreds of T Tauri and Herbig Ae star disks from obser- 
vatories such as the Spitzer Space Telescope (e.g. Furlan 
et al. 2006; Kessler-Silacci et al. 2006), the Very Large 
Telescope (e.g. van Bockel et al. 2004), the Submillimeter 
Array (e.g. Andrews & Williams 2007) and the Plateau de 
Bure Interferometer (e.g. Pietu et al. 2007). The challenge 
for theoretical astrophysicists is now not only to come up 
with a plausible model of how to overcome the various prob- 
lems in our current understanding of how planets form, but 
also to ensure that new models are consistent with the ob- 
servational data obtained from these planetary birthplaces. 

One of the challenges of theories of planet formation is 
what is often called the 'meter-size barrier'. The first stages 
of planet formation are believed to involve the coagulation 
of dust particles to aggregates of ever larger size, and even- 
tually to the size of asteroids or comets (Nakagawa et al. 
1981; Weidenschilling 1980, 1984, 1997). It was shown by 
Weidenschilling (1977) that as dust particles grow by co- 
agulation, they acquire increasingly large relative velocities 
with respect to other particles in their vicinity, as well as 
a systematic inward drift velocity as originally described 
by Whipple (1972). The high relative velocities lead to de- 



structive collisions (Blum & Wurm 2008), and thereby limit 
the growth of the aggregates to some maximum size. The 
radial drift also leads to a loss of solids toward the star 
well before large bodies can even form (Brauer et al. 2008a, 
hereafter B08a). 

Dullemond & Dominik (2005) showed that if, hypo- 
thctically, coagulation can proceed without limit, i.e., if 
no fragmentation occurs and no radial drift is present, the 
growth is so efficient that within 10 5 years very few small 
(a < 1mm) opacity bearing grains remain in the disk. The 
disk becomes optically thin and the mid- to far-infrared flux 
of the disk drops to levels well below that observed in the 
majority of gas-rich circumstellar disks around pre-main 
sequence stars. By excluding other mechanisms, Dominik 
& Dullemond (2008) argue that fragmentation is the most 
promising process for maintaining a high abundance of fine- 
grained dust in evolved disks. 

However, it has never been demonstrated explicitly that 
fragmentation is efficient enough to retain the population of 
small grains at the required levels. Moreover, if we include 
the effect of radial drift, it remains to be seen whether this 
radial drift is suppressed sufficiently by keeping the grains 
small through fragmentation. It is the goal of this Letter to 
investigate this issue and to determine how efficient frag- 
mentation should be to agree with observations. 

2. Model 

The model presented here is a combination of a ID viscous 
gas disk evolution code and a dust evolution code treating 
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the turbulent mixing, radial drift, coagulation and fragmen- 
tation of the dust. A detailed description of the model will 
follow in a later paper (Birnstiel et al., in prep.). Although 
the evolution of dust can influence the gas disk evolution, 
we neglect this effect in our model. 

The gas disk model is similar to that described in Hueso 
& Guillot (2005). The mid-plane temperature of the disk is 
approximated by following Nakamoto & Nakagawa (1994), 
taking irradiation by the central star and viscous heating 
into account. The viscous evolution of the disk surface den- 
sity I] g is described by (see Lynden-Bcll & Pringle 1974), 



dt rdr Kaa ' 
where the gas radial velocity u g is given by 

is g is the gas turbulent viscosity, 



u„ = a c. 



(1) 



(2) 



(3) 



c s denotes the sound speed, H p the pressure scale height, 
and a is the turbulence parameter (Shakura & Sunyaev 
1973). 

Grains are affected by a systematic radial inward drift 
related to headwind caused by the pressure-supported gas 
(Weidenschilling 1977), 
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where St is the Stokes number of the particle (a dimen- 
sionlcss representation of the particle size), P g is the gas 
pressure, p g is the gas volume density, and Slk is the Kepler 
frequency. In the Epstein regime, St is given by 
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where p s is the solid density of the dust grains. 

The second contribution to radial dust velocities is the 
gas drag due to the radial motion of the gas, 
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(6) 



where Sc = 1+St 2 is the Schmidt number, following Youdin 
& Lithwick (2007). 

The coupling to the gas turbulence leads to turbulent 
mixing of each dust species with a diffusion constant that 
is taken to be 

n< = 1- m 

The dust grain number density n(m, r, z) is a function of 
mass m, radius r, height above the mid-plane z, and time. 
If we define 



a(m, r) 



n(m, r, z) m 2 dz, 



(8) 



the vertically integrated time-evolution of this distribution 
can now be described by a general two-body process and 



an advection-diffusion equation, 
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Here, the right-hand side terms (from top to bottom) cor- 
respond to coagulation/fragmentation, advection, and tur- 
bulent mixing. 

Since the detailed version of the coagula- 
tion/fragmentation equation is lengthy (Weidenschilling 
1980; Nakagawa et al. 1981; Dullemond & Dominik 
2005; B08a), we include here only the general integral 
of a two-body process (first term on the right hand 
side). For a detailed description of the physics of coag- 
ulation/fragmentation, we refer to B08a; the numerical 
implementation will be discussed in Birnstiel et al. (in 
prep.). 

The physical effects that produce the relative particle 
velocities considered here are Brownian motion, relative 
radial motion, vertical settling (sec B08a), and turbulent 
motion (see Ormel & Cuzzi 2007). 

If particles collide with a relative velocity higher than 
the critical velocity Uf, they either fragment into a power- 
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law size distribution of fragments (i.e., n(m) oc m 
B08a) or the smaller body excavates mass from the larger 
one (cratering). We assume that the amount of excavated 
mass equals the mass of the smaller body. The fragmenta- 
tion velocity U{ is a free parameter in our model, and unless 
otherwise noted is taken to be 1 m/s, as suggested by ex- 
periments (e.g., Blum & Muench 1993) and by theoretical 
modeling (e.g., Stewart & Lcinhardt 2009). 

We note that the velocities and therefore also the rel- 
ative velocities are dependent on the particle size. If we 
assume that relative particle velocities are dominated by 
turbulent motion, 
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then we can estimate the Stokes number of the largest par- 
ticles by equating the relative velocity and fragmentation 
velocity, 



a ci 



(11) 



With Eq. 5, we can relate this maximum 'dimcnsionless 
size' St max to a particle size, 
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It should be noted that Eq. 10 is an approximation of the 
equations derived in Ormel & Cuzzi (2007). The detailed 
relative velocities of particles depend on the Stokes numbers 
of both particles. Hence, Eq. 12 is an order of magnitude 
estimate. 

The combined equations for coagulation, fragmenta- 
tion, and radial motion due to drift, drag, and mixing arc 
solved using the technique of implicit integration (similar to 
B08a). We use the SixPACK F-90 linear algebra package 1 
to solve the sparse matrix equation. 



available from www.engineers.auckland.ac.nz/~snor007 
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Fig. 1. Snapshots of the vertically integrated dust den- 
sity distribution of model 1 (only coagulation), model 2 
(fragmentation at 10 m/s) and model 3 (fragmentation at 
1 m/s). The grain size is given by a = (Sm/ATrps) 1 / 3 . The 
dashed line shows the evaporation radius within which no 
coagulation is calculated. The solid line shows the particle 
size corresponding to a Stokes number of unity. 



3. Results 

The simulation results presented here begin with a 0.07 Mq 
disk of power-law surface density £(r) <x r -1 (from 0.05 to 
150 AU) around a 0.5 Mq star. Unless otherwise noted, 
we used a turbulence parameter of 10~ 2 . The three mod- 
els differ in the prescription of fragmentation. In model 1, 
we accounted only for growth without fragmentation. In 
model 2, particles fragment at ut = 10 m/s, whereas in 
model 3, the fragmentation speed is taken to be 1 m/s. 
Selected results of all three models are depicted in Figure 
1. 

The top row shows the evolution in the surface density 
distribution of dust because of coagulation, radial mixing, 
and radial drift in an evolving protoplanetary disk, corre- 
sponding to model 1. The dashed line denotes the grain 
evaporation radius, which moves inwards as the tempera- 
ture falls (lower surface density caused by accretion pro- 
duces less viscous heating). The evaporation radius in our 
model is defined as the radius where the temperature rises 
above 1500 K assuming that all particles evaporate at this 
temperature. The solid line denotes the dust grain size, 
which corresponds to a Stokes number of unity. It can be 
seen that particles grow to decimeter size and then quickly 
drift inside the evaporation radius. 

The results of simulations that also include grain frag- 
mentation arc plotted in the second and third row of Figure 
1. In the innermost regions of model 2, relative particle 
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Fig. 2. Dust-to-gas ratio (solid line, linear scale) and 0.55 
fim optical depth at 10 AU (dashed line, logarithmic scale) 
as function of time for model 1 (top), model 2 (center) and 
model 3 (bottom). 



velocities are high enough (as a result of shorter dynami- 
cal times) to fragment particles already at relatively small 
Stokes number. However, the maximum size of the grain 
distribution is still strongly affected by radial drift, espe- 
cially at larger radii. Larger particles are subject to radial 
drift and therefore drift to smaller radii, where they are 
pulverized. 

In contrast, fragmentation becomes the main limiting 
factor for particle growth throughout the disk if the critical 
fragmentation velocity is 1 m/s (model 3). The result is 
that particles fragment at a small Stokes number. Particles 
with a Stokes number less than about 10 -3 arc still strongly 
coupled to the gas and not subject to strong radial drift (see 
Braucr et al. 2007). Therefore, significant amounts of dust 
can be retained for several million years if fragmentation is 
included in the calculations. 

Figure 2 compares the evolution in the dust-to-gas ratio 
(solid line) and the 0.55 pm optical depth at 10 AU (dashed 
line) of all three simulations. The optical depth is defined 
as 



cr(ra) K l/ (m)dm. 



(13) 



It can be seen that the dust-to-gas ratio declines dramat- 
ically in the first two cases. The optical depth in model 1 
drops on even shorter timescales than the dust-to-gas ratio 
because not only the dust loss but also the dust growth 
causes the optical depth to decrease. 

model 2 can retain the optical depth longer because 
small particles remain (due to fragmentation), providing 
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enough surface to keep the disk optically thick for about 0.9 
Myr. The initial dip in the optical depth is caused by the ini- 
tial condition: since initially, all mass is in small grains, the 
optical depth decreases as particles grow to larger sizes un- 
til fragmentation sets in, which increases the optical depth 
until a semi-equilibrium of growth and fragmentation is 
reached. 

If particles are already fragmented at 1 m/s (model 3), 
the dust-to-gas ratio declines only slightly, after 2 My, only 
about 4% of the dust being lost (compared to 96% and 94% 
in models 1 and 2). The optical depth decreases on longer 
timescales than in model 2 since the dust mass is lost only 
on the viscous timescale. 

This mechanism of grain retention relies on the particles 
remaining small due to fragmentation. From Eq. 11 and 12, 
it follows that the maximum size or Stokes number depends 
strongly on the critical velocity Uf and less strongly on a. If 
Uf is 10 m/s instead of 1 m/s, the corresponding particle size 
is 100 times larger meaning that the largest particles are 
now also experiencing the radial drift barrier and therefore 
drifting quickly towards the central star. 

4. Discussion and conclusion 

We have shown that grain fragmentation plays an impor- 
tant role, not only in keeping the disk optically thick (by re- 
plenishing small grains) but also in maintaining a relatively 
high abundance of solids of all sizes. If particles fragment at 
relatively small sizes, they can stay below the size-regime 
where they are affected by radial drift. 

The simulation results presented here suggest that frag- 
mentation of relatively low critical velocity (a few m/s) 
is needed to retain the dust required for planet forma- 
tion, although the formation of planetesimals probably re- 
quires periods (and/or regions) of quiescence such as dead 
zones or pressure bumps (see Kretke & Lin 2007; Brauer 
et al. 2008b) . These locally confined quiescent environments 
could decrease both the relative and radial velocities of 
particles, thus allowing particles possibly to either grow 
to larger sizes or to accumulate until gravitational insta- 
bilities become important (see Johansen et al. 2007; Lyra 



et al. 2009). These effects are not yet included in this model 
and will be the subject of future research. 

Assuming that the proposed mechanism of dust reten- 
tion is effective in protoplanetary disks, we can constrain 
the critical fragmentation velocity to be smaller than 10 
m/s since velocities > 10 m/s cannot explain the observed 
lifetimes of disks even in the case of small a. 
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